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We study the spectra of gamma ray signals that arise from dark matter annihilation in the 
universe. We focus on the large class of theories where the photon spectrum includes both continuum 
spectrum of gamma rays that arise from annihilation into Standard Model states at tree level, as 
well as monochromatic gamma rays arising from annihilation directly into two photons at the one 
loop level. In this class of theories we obtain lower bounds on the ratio of the strength of the gamma 
ray line relative to the gamma ray continuum as a function of the dark matter mass and spin. These 
limits are obtained from the unitarity relation between the tree level amplitude of the primary 
annihilation channel and the imaginary part of the loop level amplitude for annihilation directly 
into photons, with the primary decay products running in the loop. These results are exact in the 
limit that dark matter annihilation is exclusively to a single Standard Model species, occurs through 
the lowest partial wave and respects CP. Away from this limit the bounds are approximate. Our 
conclusions agree with the known results in the literature in the case of the Minimal Supersymmetric 
Standard Model (MSSM). We use the Fermi-LAT observations to translate these limits into upper 
bounds on the dark matter annihilation cross section into any specific Standard Model state. 



I. INTRODUCTION 

Cosmological measurements have now established that 
about 80% of the matter in the universe is composed of 
non-luminous, non-baryonic dark matter [1]. However, 
the precise nature of the particles of which dark matter 
is composed remains a mystery. Weakly Interacting 
Massive Particles (WIMPs) - thermal relics with weak 
scale masses and weak scale cross sections with visible 
matter - constitute one well-motivated class of dark 
matter candidates. In these theories the relic abundance 
of the WIMP is set by its annihilation cross section into 
Standard Model (SM) particles, and turns out to be 
naturally of the right order to explain observations. 

One approach to dark matter detection involves search- 
ing for the products of WIMP annihilation in the uni- 
verse, such as photons. Dark matter is constrained 
to be neutral under electromagnetism, and therefore 
in renormalizable theories WIMPs cannot annihilate 
directly into photons at tree level. Nevertheless, a 
continuum spectrum of photons arises from decays of 
the primary annihilation products, and also from final 
state radiation off charged final states. The spectrum 
also includes monochromatic gamma rays that arise from 
annihilation directly into 77, 7Z and jH final states, 
usually at loop level. 

The continuum spectrum of photons arising from dark 
matter annihilation into any specific final state is to 
a large extent independent of the detailed form of the 
matrix element for the process, depending only on the SM 
quantum numbers of the particles in the final state and 
the dark matter mass. This allows limits on the observed 
diffuse gamma ray flux from regions of high dark matter 
density to be translated into robust bounds on the rate 



of dark matter annihilation into any such state. 

Over the last few years the Large Area Telescope 
(LAT) aboard the Fermi Gamma-ray Space Telescope 
has been making precise observations of the gamma ray 
sky. Limits on the rate of dark matter annihilation 
into various final states have been placed by Fermi- 
LAT measurements of the gamma ray flux from nearby 
galaxy clusters [2], dwarf galaxies [3, 4], the Galactic 
center [5] and subhalos [6]. The diffuse gamma ray flux 
measurement from portions of the sky [7, 8] as well as 
the isotropic near full-sky [9, 10] have been used to set 
bounds on dark matter annihilation. 

In typical models of WIMP dark matter, the relic 
abundance is set by tree level annihilation to final 
states consisting of two SM fermions or two weak 
gauge bosons. In this scenario monochromatic gamma 
rays corresponding to annihilation to the 77, jZ and 
-fH final states are only generated at one loop. The 
strengths of the corresponding gamma ray lines are 
highly model dependent, being very sensitive to the 
detailed structure of the corresponding matrix element, 
which in turn depends on the details of the vertices in 
the theory. Furthermore, the amplitude will in general 
receive significant contributions from unknown non-SM 
states (associated with new physics) running in the 
loop. Therefore, the Fermi-LAT bounds on gamma ray 
lines [11] have primarily been used to probe specific well- 
motivated models. In particular, WIMP models of this 
type for which the one loop cross sections to the 77 
and 7Z final states have been computed include the 
Minimal Supersymmetric Standard Model (MSSM) [12- 
15], Universal Extra Dimensions (UED) [16, 17], the 
Littlest Higgs Model with T-Parity [18, 19] and the Inert 
Doublet Model [20]. In theories where dark matter is 



a scalar or Majorana fermion, annihilation to the 
final state is suppressed in the non-relativistic limit. 
However, the cross section to this and to the 7Z final 
state have been computed in a specific theory of Dirac 
dark matter [21] that arises in a class of Randall- 
Sundnim models. A model-independent approach based 
on effective operators that correlates the signal in direct 
detection experiments with the approximate line strength 
has also been developed [22], (see also [23]). 

In this paper, we place lower bounds on the strengths 
of gamma ray lines from dark matter annihilation that 
are applicable to a large class of theories. We focus 
on models where the photon spectrum includes both 
continuum gamma rays that arise from annihilation 
into two body SM final states at tree level, as well as 
monochromatic gamma rays arising from annihilation 
directly into two photons at the one loop level. For this 
class of theories we obtain lower limits on the ratio of the 
strength of the gamma ray line relative to the gamma 
ray continuum as a function of the dark matter mass 
and spin. These limits are obtained from the unitarity 
relation between the tree level amplitude of the primary 
annihilation channel and the imaginary part of the loop 
level amplitude for annihilation directly into photons, 
with the primary decay products running in the loop. 
These bounds are exact in the limit that dark matter 
annihilates exclusively into a single SM species, that the 
annihilation is dominated by the i = partial wave and 
respects CP. Relaxing these assumptions renders our 
limits approximate, but does not invalidate them. In 
particular, we stress that these results do not depend on 
the specific form of the coupling between the dark matter 
particle and the SM field running in the loop, and are 
in that sense model-independent. Our conclusions agree 
with the known results in the literature in the case of the 
Minimal Supersymmetric Standard Model (MSSM). We 
use the Fcrmi-LAT observations to translate these limits 
into upper bounds on the dark matter annihilation cross 
section into any specific Standard Model state. 



where the sum over \X) runs over all final states into 
which the dark matter particle can annihilate. 

In the center of mass frame, two particle states can 
be labelled by the internal quantum numbers of the 
particles, the total energy, and the angular momentum 
quantum numbers \J,M]L,S), where J is the total 
angular momentum, M is the component of total angular 
momentum along any fixed axis, L is the orbital angular 
momentum and S the total spin angular momentum. 
Dark matter annihilation in haloes occurs in the highly 
non-relativistic regime. We therefore focus on the case 
where annihilation occurs through the L = partial 
wave, so that the initial state \i) has J = S. While J 
and M are conserved in the annihilation process, this is 
not true of L or 5. 

There is an alternative basis for labelling two particle 
states in the center of mass frame, the helicity basis, 
which will also prove useful. In this basis, in addition 
to the internal quantum numbers of the particles and 
their total energy, the states are labelled by the angular 
momentum quantum numbers |J, M;Ai,A2), where Ai 
and A2 are the helicities of the two particles. Unlike 
J and M, the helicities are not conserved in the an- 
nihilation process. The transformation that relates the 
I J, M; Ai, A2) basis to the | J, M; L, S) basis may be found 
in the classic paper by Jacob and Wick [24]. 

Now, if the theory is invariant under time reversal 
(or equivalently CP), and the states |z) and [/) are 
eigenstates of angular momentum, labelled either by 
\J,M;L,S) or by |J, M;Ai,A2), then, as shown in the 
appendix, (/|T|i) = {i\T\f). Eq. (3) then simplifies to 



2Im(/|T|i) = ^(/|Tt|X)(X|T|i). 



(4) 



X 



At lowest order in perturbation theory, this equation 
relates the imaginary part of the loop amplitude for 
annihilation into the two photon final state to the tree 
amplitude for annihilation into two body SM states. 
Squaring, we get 



II. LOWER LIMITS ON LINE STRENGTHS 

In this section, we explain in detail how these limits are 
obtained. Unitarity implies that the 5-matrix satisfies 



S^S = 1. 



(1) 



Writing S = 1 + iT, where T is the transition matrix, we 
obtain 



i(T - Tt) = T^T. 



(2) 



Consider the matrix element of this equation between the 

initial state \i) consisting of the two dark matter particles, 
and the final state [/) consisting of two photons, 

-i{f\{T-T^)\i) = ^{f\T^\X){X\T\i), (3) 

X 



4|Im(/|T|i)|^ 



X:(/|Tt|X)(X|T|z) 



X 



(5) 



We now restrict ourselves to the case where dark 
matter annihilation at tree level is exclusively to a single 
SM species. If the SM state \X) is further characterized 
by fixed values of L and S (or fixed values of Ai and A2) 
for any given values J and M of the initial state then 



4|im(/|r|z)|2^|(/|rt|x)|2|(x|r|z)|^ 



(6) 



This equation implies that for all initial states \i) that 
annihilate at tree level exclusively to the SM state \X) 
characterized by fixed L and S (or fixed Ai and A2), the 
ratio 



|Im(/|T|z) 



\\{f\T\X)\' 



(7) 
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FIG. 1. Model independence of the lower bound on the ratio of loop and tree annihilation cross sections shown graphically. 
The X particle is a standard model state, and x is the dark matter particle. For purposes of obtaining a limit, dark matter 
annihilation can be modeled by the decay of a boson, 



is a constant that depends only on SM parameters, and is 
otherwise independent of Furthermore, under these 
conditions the quantity 



Ml^>\lf\T\X)\^ 



(8) 



is bounded from below. The numerator of the expression 
on the left hand side is proportional to the cross section 
for \i) — > I/), while the denominator is proportional to 
the cross section for \i) — >• \X). The right hand side is 
proportional to the absolute value of the tree level matrix 
element for annihilation of the state \X), consisting 
of two SM particles, into two photons. This can be 
calculated within the SM. Therefore this expression can 
be translated into a lower bound on the rate for xx ^ 77 
relative to the rate for xx X X. To obtain the bound 
we proceed as follows. 

• We construct a theory with a bosonic particle $ of 
mass m$ = 2my^ that decays exclusively into the 
state |X) at tree level, and into two photons at one 
loop. The spin and CP properties of $ are chosen 
to be identical to those of the initial state composed 
of two dark matter particles. 

• We calculate the contribution to the decay rate 

for $ — > 77 that arises purely from the imaginary 
part of the corresponding one loop amplitude. We 
denote this quantity by rim($ 77). 



• We calculate the decay rate for $ 

by r($ ^ X X). 



X X, denoted 



As shown above, the ratio of these two decay rates 
will be independent of the details of interactions 
of Therefore, the decays of <& can be used to 
model dark matter annihilation. Hence, we obtain 
the bound (see Fig. 1): 



o-(xx ^ 77) ^ rim($ -j> 77) 



a{xX^XX) - T{^^XXy 



(9) 



• Using the fact that the continuum spectrum of 

gamma rays corresponding to dark matter annihi- 
lation to the state \X) is known, this result can be 



translated into a lower bound on the strength of 
the gamma ray line relative to the continuum. 

For the bound to apply, dark matter annihilation at 
tree level must occur exclusively into a single SM final 
state. Furthermore, this state must be characterized 
by definite values of L and S (or definite Ai and A2), 
for any given J and M . Under what circumstances are 
these criteria satisfied? This turns out to depend on the 
dark matter spin, on whether the dark matter particle 
is its own anti-particle, and on the masses ands spins of 
the SM particles in the state \X). We now consider the 
various possibilities in turn. 

Scalar dark matter 

Consider first the case where dark matter is a scalar, 
either real or complex. Then the initial state is CP 
even, and has J = 0. Angular momentum conservation 
implies that annihilation into light fermions is helicity 
suppressed. Therefore we expect that annihilation will 
occur primarily to the heaviest SM ferniion species that 
is kinematically accessible, or alternatively to W~^W~ 
(excluding neutral states which obviously do not con- 
tribute to photon line signal at one loop). In the case 
of annihilation to fermions, angular momentum and CP 
conservation imply that the final state must have L = 
1,S = 1, while L = 0, 5 = is forbidden. Therefore a 
bound can be obtained for this annihilation mode. 

For annihilation to W^W~ on the other hand, 
the conservation laws allow both i = 0, 5 = and 
L = 2, S = 2 final states, forbidding only L = 1,5= 1. 
Therefore our formalism is not directly applicable. 
Nevertheless, in the next section we shall see that in 
various kinematic limits a bound can indeed be obtained. 

Majorana fermion dark matter 

We move on to the case where dark matter is 
a Majorana fermion. Since Majorana fermions are 
identical particles, anti-symmetry of their wave function 
implies that if L = 0, S is also then zero so that J = 0. 
However, unlike the case of scalar dark matter, the 
initial state is now CP odd. As before annihilation 
to light fermions is disfavored by angular momentum 



3 



considerations, and so the heavy fermion and W~^W~ 
final states are again expected to dominate. Consider 
first annihilation to fermions. Angular momentum and 
CP conservation imply that the final state must have 
L = 1,S=1, while L = 0, 5' = is forbidden. Therefore 
our formalism applies to this annihilation channel. What 
about the W~^W~ channel? Now the conservation laws 
allow only the L = 1, S = 1 final state, while forbidding 
the X = 0, 5 = and L = 2, S = 2 final states. Therefore 
our formalism applies to this channel as well. 



Dirac fermion dark matter 

If dark matter is a Dirac fermion there are two 
possibilities for the total angular momentum of the 
initial state, J = or J = 1. Annihilation will in general 
proceed through both these channels, and since J is a 
conserved quantum number there is no interference. The 
bound in each channel can be calculated independently. 
The weaker of these two limits is then the true bound, 
corresponding to the case where annihilation occurs 
entirely through that channel. For J = 0, CP is odd 
and the analysis is identical to that of Majorana fermion 
dark matter considered above. However, if J = 1, 
annihilation to the 77 final state is forbidden by the 
Landau- Yang theorem [25, 26]. Therefore dark matter 
can annihilate to SM final states through this channel 
without giving rise to a line signal from the two photon 
final state. We conclude that in the case that dark 
matter is a Dirac fermion, no general bound is possible. 



III. COMPUTATION OF LIMITS 

In this section we place lower limits on line strengths 
by explicitly calculating the ratio of the decay rates of 
the boson $: 



Tim ('I' 11) 

T{^^x xy 



(10) 



As outlined earlier, different choices of the spin and CP 

properties of the boson <1> map on to different dark matter 
candidates. We choose the mass of $ to be 2m^, to 
ensure that the intermediate state particles have exactly 
the same total energy as in dark matter annihilation. 
This is required so that the quantum numbers of the 
intermediate state can be exactly the same in the two 
cases, which is required for the mapping to go through. 
The final result will in general be seen to depend on the 
velocity /3 of the intermediate state particle X, 



/3 = 1 1-^ 



Scalar Dark Matter 



(11) 



In the case where the dark matter particle is a scalar, 

the incoming state is restricted to be in a J = state 
for non-relativistic annihilation. The initial state is then 
a CP even state with zero total angular momentum, 
which is duplicated when the boson $ is a scalar particle, 
which we label by (j). 



Vector boson dark matter 

Finally we consider the case where dark matter is a 
real vector boson. Since the initial state is composed of 
identical particles, the wave function must be symmetric, 
allowing J = or J = 2, but forbidding J = 1. Both 
J — Q and J — 2 are CP even. For J — the analysis 
is identical to that of scalar dark matter considered 
above. However, for J = 2 a separate analysis is needed. 
Annihilation may occur to either light or heavy fermions, 
or to TF+T'T'"^. In particular, annihilation to light 
fermions is no longer disfavored by angular momentum 
conservation. Annihilation to fermions may occur either 
through any of the L = {1,2, 3}, 5 = 1 channels. 
Therefore, even if annihilation occurs exclusively to a 
single fermion species, it is not possible to obtain a 
general bound. However, it is possible to obtain a bound 
in various kinematic limits. 

Annihilation to W~^W~ can also occur through mul- 
tiple channels, including L = 2,S = 0, and L = 
{0, 1, 2, 3, 4}, S* = 2. Therefore, in general it is not 
possible to constrain this channel either. However, we 
shall see that it is possible to obtain a bound in the limit 
that the dark matter mass is close to the W mass, so that 
the final state W bosons are non-relativistic. 



Scalar dark matter annihilation to fermions 

As explained in the introduction, symmetry consider- 
ations require the fermions to be in the L — I, S — 1 
final state. We assume that dark matter annihilates 
predominantly to a single species of SM fermions //. 
In the absence of large new sources of chiral symmetry 
breaking in the theory, this is a safe assumption, since the 
matrix element for annihilation a given fermion species is 
proportional to the fermion mass. Then annihilation to 
heavier fermions such as tops and bottoms are preferred, 
provided those channels are kinematically open. 

We begin from the following interaction Lagrangian, 
which is assumed to be valid below the scale of elec- 
troweak symmetry breaking: 



(12) 



Prom this Lagrangian it is straightforward to calculate 
the ratio of the imaginary part of the decay to two 
photons to the total decay rate to //. We obtain for 
the bound 



77) NcQ^e'^mj 



/3 [tanh-^/3]^ (13) 



where Q is the electric charge of the fermion (e.g. Q = § 
for top quarks) and A^c = 3 is the color factor to be 
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Dark Matter 


Initial spin 


Annihilation 


Bound 


Channel 


Mode 


Scalar 


J = 


WW 


L = 0,5 = 
L = 2,5 = 2 


In NR / UR limits. 






ff 


L = 1,5 = 1 


/ 


Majorana Fermion 


J = 


WW 


L = 1, 5 = 1 


/ 




ff 


L = 0,5 = 


/ 




J = 


WW 


L = 1,5 = 1 


/ 


Dirac Fermion 




ff 


L = 0,5 = 


/ 




r — 1 


Forbidden 




J = 


WW 


L = 0,5 = 
L = 2,5 = 2 


In NR / UR Umits. 


Real Vector Boson 




ff 


L = 0,5 = 


/ 




J = 2 


WW 


L = 2,5 = 
L={0,1,2,3,4},5 = 2 


In NR limit. 




ff 


L = {1,2,3},5 = 1 


In NR / UR Umits. 



TABLE I. Summary of the bounds corresponding to the different dark matter candidates and annihilation channels. Wc have 
also indicated cases where a bound is only applicable when the tree-level annihilation products (on-shell intermediate states in 
the loop) are non-relativistic (NR) or ultra-relativistic (UR). 



included when the fermionic states are quarks. 

Scalar dark matter annihilation to W bosons 

As explained carlicir. in this case annihilation can 
proceed through either the L = Q,S = Q channel or 
through the L = 2,S = 2 channel. Therefore it is only 
possible to obtain model independent bounds in specific 
kinematic limits. 

Consider first the limit that the dark matter mass 
is close to the W mass, — mw ^ m-w- Then 
the annihilation products will be non-relativistic and we 
expect that the L = 0,S = final state will dominate. 
Therefore, in this limit a bound can be obtained. As 
before we model the annihilation by the decay of a scalar 
particle (p. We choose the coupling of to bosons to 
have the simple form 



(14) 



and calculate the decay rates to WW at tree level, and 
to 77 at loop level. The required ratio of decay rates is 
given by 



rim(<^ 77) 



3e4 



r/3 



T{(t)^WW) 647r2' 

(non-relativistic limit). 



(15) 



This is the bound in the non-relativistic limit. 

In the opposite limit where the dark matter mass is 

much larger than the W mass, ^ ttiw, the W 
bosons are ultra-relativistic. Then the Goldstone boson 



equivalence theorem applies, and the longitudinal and 

transverse components of the W bosons correspond to 
distinct physical states. Annihilation can occur either 
to two longitudinal W bosons or to two transverse W 
bosons, and it is reasonable to assume that there is no 
large cancellation in the contributions to the 77 ampli- 
tude from these different states. Wc therefore consider 
them separately. (Annihilation to one transverse and 
one longitudinal W is forbidden by angular momentum 
conservation, and so need not be considered.) 

If annihilation is exclusively to longitudinal W^'s, then 
there is a unique final state labelled by |0, 0,0,0) in 
the helicity basis |J, M, Ai,A2), and our formalism can 
immediately be applied. If annihilation is exclusively to 
transverse Ws, then angular momentum conservation 
allows both the |0, 0, -h, -h) and |0,0, — ,— ) final states. 
However, only one linear combination of these states is 
CP even, while the other is CP odd. Since annihilation 
can only occur to the CP even state, once again our 
formalism is applicable. 

We first consider annihilation to longitudinal W 
bosons. We couple the scalar to the Higgs doublet, 
and consider decays into charged Higgses, and into 
two photons through a loop of charged Higgses. The 
Goldstone boson equivalence theorem guarantees that in 
the ultra-relativistic limit, the rates for these processes 
are identical to the rates for the corresponding processes 
involving W bosons. The interaction Lagrangian takes 
the form 



Cint = a(f> H^H. 



(16) 
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The ratio of the decay rates in this case turns out to be 
suppressed. It scales roughly as the following, 



rim(0 77) 
r{(j)^ww) 



m 



w 



log 



Ami 



(17) 



(ultra-relativistic limit). 



Since this vanishes in the limit 



m^/m^ 0, this shows 
that there is in fact no bound in the case of annihilation 
to longitudinal W bosons in the ultra-relativistic limit. 

We now consider annihilation to the transverse polar- 
izations. We start from the interaction Lagrangian 



(18) 



In the ultra-relativistic limit, this interaction leads to 
decays of the scalar (jj exclusively into transverse polar- 
izations of the W boson. The ratio of cross-sections is 
obtained as 



rim(0 77) 



WW) 327r2 



log 



4to^ 
ml, 



-\ 2 



(19) 



(ultra-relativistic limit) . 



This is then the lower bound in the case of annihilation 
exclusively into transverse polarizations of the W. 

In general, annihilation will occur into both transverse 
and longitudinal polarizations. From the discussion 
above, it follows that if Ft is the branching fraction into 
transverse polarizations of the W in the ultra-relativistic 
limit, the bound is given by 



rim(0 77) 
r{(j) WW) 



loe 



4m: 



327r2 

(ultra-relativistic limit) 



(20) 



This formula is valid in the ultra-relativistic limit pro- 
vided 



Ft > 



mt 



(21) 



If the branching fraction to transverse polarizations is 
smaller than this, then there is no bound in the limit 
m^/m^ 0. 



B. Majorana Fermion Dark Matter 

In the case when the dark matter particles are 
Majorana fermions, the initial state again has total 
angular momentum J — 0. This is because in the 
non-relativistic limit L — is picked out, and then 
5 = is required by the overall antisymmetry of the 
total wavefunction of the two fermion state. This 
configuration is CP odd. The angular momentum and 
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FIG. 2. Spectral flux E^d^/dE in the continuum relative to 
the lower limit for the line obtained from our calculation. 
We show results for annihilation into WW (upper panel) 
and bb (lower panel) for a Majorana fermion dark matter. 
The annihilation cross section was taken to be {aAv) 



3 X 10" 



cm^ s"^ and = 100 GeV. The width of the line 



is set by the Fermi-LAT resolution at 100 GeV, 11% FWHM. 



CP quantum numbers of the initial state are exactly 
those of a pseudo-scalar particle. Therefore we can 
calculate the ratio by decaying a pseudo-scalar particle, 
which we denote by tp. 

Majorana dark matter annihilation to fermions 

We begin by considering annihilation to fermions. 
Once again we assume annihilation is exclusively to a 
single SM fermion species. As in the case of scalar dark 
matter this is a safe assumption, since annihilation to 
light fermions is chirality suppressed in the absence of 
large new sources of chiral symmetry breaking in the 
theory. 

Angular momentum conservation and CP symmetry 
require the final state fermions to be in the L = 0, S = 
configuration rather than L = 1,5 = 1. Our 
starting point is the following interaction Lagrangian 
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which couples the pseudo-scalar y> to the fermions //: 
jCint=iXfl^f<P- (22) 

Prom explicit calculation we obtain the bound in this case 
as 

where as before Q denotes the electric charge, and Nc 
denotes the color factor associated with the fermion. 
We have verified that this agrees with the result in 
the literature for the case of neutralino annihilation 
to fermions in the Minimal Supersymmetric Standard 
Model (MSSM). 

Majorana dark matter annihilation to W bosons 

We move on to considering annihilation to W bosons. 
The only angular momentum quantum numbers for 
the WW final state consistent with J = and CP 
conservation are L = 1, S' = 1. As before, wc model this 
annihilation by pseudo-scalar decay. The pseudo-scalar 
if is coupled to the SM in the following way: 



Ant = J^Tl:{F^,F^"'). 



This leads to the bound: 



(24) 



(25) 



This agrees with the result in the literature for neutralino 
annihilation to W bosons in the MSSM. 

Fig. 2 shows the lower limits on the strengths of the 
gamma ray line relative to the gamma ray continuum if 
Majorana fermion dark matter annihilates exclusively to 
the bb and WW final states respectively. The widths of 
the emission lines are drawn so as to correspond with the 
resolution of the Fermi-LAT, at 11% of FWHM. 



C. Real vector boson dark matter 

Real vector boson dark matter annihilation can pro- 
ceed through either the J = channel or the J = 2 
channel. The J = 1 channel is not allowed because 
the wavefunction is required to be symmetric under 
interchange of the two identical dark matter particles. 
Both the J = and J = 2 initial states are CP even. 

If dark matter annihilation occurs primarily through 
the J = channel, then the bounds from scalar dark 
matter annihilation to the corresponding SM states apply 
here as well. On the other hand, if the primary mode of 
annihilation is through the J = 2 channel, then obtaining 
the bound involves calculating the decays of a spin-2 
CP even particle. The model independent lower limit 
on the strength of the line corresponds to the weaker of 
the bounds obtained in the J = and J = 2 cases. 



Wc use the decays of a massive spin-2 graviton to 
model the J = 2 annihilation process. The massive 
graviton is taken to couple to the stress energy tensor 
of the matter fields, in analogy with the couplings of 
the Kaluza-Klein graviton in extra-dimensional theories 
[27, 28]. 

Vector (lark mat tor uuuiliilaiiou to fermions 

We begin by considering annihilation to fermions. The 
bounds we obtain assume annihilation exclusively to a 
single species of SM fermion. Away from this limit our 
bounds are approximate. The two fermion final states 
that can arise from annihilation of the J = 2 initial state, 
after taking into account angular momentum and CP 
conservation, have angtilar momentum quantum numbers 
L = 1,5 = 1, L = 2,5 = 1 and L = 3,5 = 1. 
This multiplicity of available states means that it is only 
possible to obtain a bound in specific kinematic limits. 

We first consider the limit where the dark matter mass 
is close to the mass of the fermion species it annihilates 
into, so that the outgoing fermions are non-relativistic. 
Then the L = 1, S = 1 final state dominates, and our 
formalism is applicable. 

Consider a massive Dirac fermion / that couples to a 
massive graviton h^n, . The term in the Lagrangian which 
is relevant for the on-shell graviton decay is the following: 



(26) 



The coupling constant k has dimensions of inverse mass. 
The mass of the fermion couples to the trace of the 
graviton /iJJ, which vanishes for the on-shell decay. The 
limit in this case is obtained very weak because of the 
p-wave suppression of the amplitude: 



rim(/l 77) 



r(/i ^ //) 



J=2 



1207r2 

(non-relativistic limit) . 



(27) 



This is significantly weaker than the corresponding 
bound on decays from the J = state. We conclude 
that for real vector dark matter annihilation to heavy 
fermions in the non-relativistic limit, it is the bound from 
the J = 2 initial state that applies. 

Next we consider the case of annihilation to very light 
SM fermions, such as electrons or muons. In particular, 
we work in the limit where the dark matter mass is 
much larger than the mass of the final state fermions. 
Then the final state fermions are ultra-relativistic, and 
can be treated as massless. In this limit, left- and 
right-handed fermions of the same flavor are effectively 
different SM species, and should be considered separately. 
We now show that in the ultra-relativistic limit, for 
annihilation exclusively to a single SM fermion species 
of definite chirality, it is possible to obtain a bound. To 
understand how this arises we work in the helicity basis, 
and consider annihilation exclusively to left-handed Weyl 
fermions (and right-handed anti-fermions) . The unique 
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final state is then | J, M; +i, —5), and so our formalism 
can immediately be applied. 

The couplings of a massless chiral Weyl fermion / to a 
massive graviton take the form: 



We obtain the result 
rim(/i 77) 



7=2 



1447r2 

(ultra-relativistic limit). 



(28) 



(29) 



Since annihilation to massless fermions from the initial 
state with J = is forbidden, this is the applicable 
bound. 

In this scenario, in contrast to the cases of scalar and 
Majorana fermion dark matter, the annihilation cross 
section to SM fermions is not chirality suppressed. It 
is then very natural for dark matter to couple identically 
to all three flavors of a SM (chiral) fermion species and 
annihilate with equal strength to all of them. In the 
absence of other annihilation modes, Eq. 29 gets modified 
to 



Tim (ft H> 77) 



r(/i ^ //) 



4^4 



J=2 



1447r2 

(ultra-relativistic limit) 



(30) 



where Nf is the number of fiavors. 

Vector dark matter annihilation to W bosons 

For annihilation to W bosons from the J = initial 
state, the earlier bounds from the case of scalar dark 
matter apply in the appropriate kinematic limits. In 
the case of the J = 2 initial state, even after imposing 
angular momentum and CP conservation, there are still 
several channels that can contribute. It is possible to 
obtain a bound in the limit that the dark matter mass 
is close to the W boson mass, — mw ^ mw, when 
annihilation to the L = 0, S = 2 final state dominates. 

The relevant part of the Lagrangian, in unitary gauge, 
takes the form: 



-m^^W+W-] + fx ^ u) (31) 



Explicit calculation leads to the bound: 
rim(/i^77) o 



T{h WW) 



J=2 



207r2'^ 
(non-relativistic limit). 



(32) 



This bound is (slightly) stronger than than that obtained 
for annihilation through the J = channel (Eq. 15). We 
conclude that the J = channel provides the appropriate 
bound for vector dark matter annihilation to WW in the 
non-relativistic limit. 



IV. CURRENT AND FUTURE 
EXPERIMENTAL LIMITS 

In this section we combine these results with the 
current experimental constraints on gamma ray lines 
obtained by the Fermi-LAT collaboration [11] to place 
upper bounds on the rate of dark matter annihilation 
into various SM final states, and estimate potential future 
sensitivities. 

We consider each of scalar, Majorana fermion and 
real vector boson dark matter in turn and consider 
annihilation into a few plausible SM final states. We have 
obtained a lower limit on the strength of the photon line 
given the annihilation cross section into each such final 
state. The null result in Fermi-LAT line search then puts 
an upper bound on this annihilation cross section. 

Our results are then compared to the corresponding 
bounds on dark matter annihilation obtained from the 
Fermi-LAT constraints on continuum gamma rays. This 
allows us to compare the relative strength of the bounds 
from the line and the continuum photon spectra searches. 

The flux in a line is given by 



where 



dE 



m 

dE 



(aAv) J dN 
87rm| Jo dE ' 



= 2d{E^ - m^) 



(33) 



(34) 



for the given annihilation channel with annihilation cross 
section {axv), and Jo = l/[8.5 kpc(0.3 GeV cm~^)^]. 
The normalized integral of the mass density squared over 
the observational regions for the line search is 



(35) 



J = io / p {rgai{b,i,x))cosb dx db ' 



where rgai(6,£,x) is the radial coordinate of the density 
distribution. The regions included in the line search 
are the Galactic caps at |6| > 10° and the 20° x 20° 
region in the Galactic center at |6| < 10° and \£\ < 10°. 
We adopt a minimal model of the Milky Way Galactic 
dark matter profile which minimizes the signal yet is 
consistent with constraints from Galactic dynamics and 
is consistent with profiles expected in cold dark matter 
halo formation. To do so, we take an Einasto profile 
for the dark matter density distribution with parameters 
that minimize J and are consistent with the dynamical 
constraints in Ref. [29]. The Einasto profile is 



PEinasto('") = Ps exp 



1 



(36) 



and the parameters from Ref. [29] within 68% confidence 
level (CL) that give a lower J are a = 0.22, r,, = 
21 kpc, with r0 = 8.28 kpc and rg determined from the 
local solar dark matter density Pq = 0.385 GeV cm^'^ 
The value of the normalized density squared toward the 
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FIG. 3. Upper, middle and lower rows of panels show current and future 95% CL constraints on dark matter annihilation 
channels in the scalar, Majorana, and vector dark matter cases respectively. The colored regions are from constraints from 
the isotropic DGRB on Galactic and extragalactic dark matter annihilation, as described in the text. The diagonally hatched 
regions are from current constraints from Fermi-LAT on the minimal line contribution in the respective channels, as discussed 
in the text. The vertically hatched regions are forecasts for the limits after 10 years of observation with Fermi-LAT, assuming 
Poisson statistics at the current energy resolution. In the case of scalar dark matter to WW (upper right panel), the solid line 
indicates the true limit for ultra-relativistic WW production, and the non-outlined hatched region is an interpolation of the 
limit to the non-relativistic regime at 171^ ~ 100 GeV. 



Galactic caps in this case is Jcap — 28.6 and toward 
the Galactic Center's 20° square is Jgc = 37.0. These 
limits depend on the adopted profile among dynamically 
consistent models at the level of ^-^50%, as shown in 
Ref. [11]. Our adopted profile and parameters give a 
minimal J^, and therefore provide a conservative limit, 
yet are consistent with cold dark matter halo structure 
formation. 

In the case of Majorana fermion dark matter we 
consider annihilation to the bb, t^t^ and WW final 
states. For scalar dark matter we consider exactly the 
same final states, but with the added restriction in 
the WW case that the I^-bosons are ultra-relativistic 
and transverse. In the case of real vector boson dark 
matter, we first consider the case where annihilation 
occurs equally to all three generations of ultra-relativistic 



right-handed down-type quarks 66, ss and dd, then the 
case where it occurs equally to all three generations of 
ultra-relativistic right-handed charged leptons rf, ij.fl 
and ee and finally the case where annihilation is to non- 
relativistic M^-bosons. The 95% confidence level (CL) 
limits are shown in Fig. 3. 

The continuum limits, as shown in Fig. 3, are derived 
as in Ref. [9] from Galactic and extragalactic contri- 
butions to the isotropic diffuse gamma ray background 
(DGRB) as measured by Fermi-LAT [30], with the 
exception that here we take the substructure boost to 
be a very conservative B = 2.3 within the Galactic 
halo component of the DGRB contribution. This boost 
value is found from the substructure enhancement as 
in Ref. [9], with minimal substructure parameters, and 
with a partial cancellation of the total luminosity boost 



9 



contribution from our position within the Galactic halo. 

The minimal line limits arc in all cases weaker than 
continuum limits at this time from Fcrmi-LAT's measure- 
ment of the isotropic DGRB. Note that limits from the 
stacking of dwarf galaxies are stronger than that from the 
DGRB by a factor of ~10 [4] and the isotropic DGRB 
sensitivity may be enhanced [31]. We compare to the 
isotropic DGRB constraints here because they are the 
most conservative among the annihilation constraints. In 
Fig. 3 we also show the expected future 10-year Fermi- 
LAT line limits assuming enhancement of the sensitivity 
with the exposure time (t) of Fermi-LAT as The 
line sensitivity for Fermi-LAT may improve substantially 
with systematic improvement at better than the Poisson 
count rate. For example, this may occur with enhanced 
energy resolution for lines over the lifetime of Fermi-LAT. 
Moreover, for certain cases of dark matter spins, any 
future gamma-ray experiment with even higher energy 
resolution (e.g. GAMMA-400 [32]) may prove to have 
line sensitivities or constraints on dark matter annihila- 
tion be more significant than that for the continuum. 

V. CONCLUSIONS 

In models where the dark matter particle has a primary 
annihilation channel to a unique final state in the SM, we 
have derived a robust lower limit on the cross section of a 
loop-induced annihilation mode to photons, for a variety 
of choices of spin of the dark matter particle as well as 
the SM states it annihilates into. This bound is based on 
unitarity considerations that relate the imaginary part of 
the loop induced amplitude for annihilation to photons 
to the cross section of the primary annihilation chan- 
nel. Since the spectrum of continuum photons emitted 
through bremsstrahlung from the primary annihilation 
products can be reliably computed, this lower bound also 
relates the minimal strength of a gamma ray line from 
dark matter annihilation to the size of the continuum 
spectrum. While the very conservative bounds for a 
gamma ray line obtained in this way are less stringent 
than the bounds from the continuum photons, these 
results help identify combinations of initial and final 
state quantum numbers for which a model-dependent 
calculation of the full amplitude (rather than only its 
imaginary part) for annihilation into photons can easily 
give significantly stronger bounds than those derived 
from continuum emission. 
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Appendix A: Time-Reversal Invariance and the 
T-Matrix 

In this appendix we show that if the theory respects 
time-reversal invariance, then the T-matrix is symmetric 
in the angular momentum basis, 

{f\T\i) = {i\T\f) , (Al) 

where \i) and |/) are angular momentum eigenstates. 

The fact that the theory is time-reversal invariant 
implies that 

{f\T\i) = {ei\T\Qf) , (A2) 

where O is the time-reversal operator. The action of 
time-reversal on eigenstates of angular momentum is 
given by [24] 

Q\J,M;L,S) = {-iy-'^\J,-M;L,S) 
e|J,M;Ai,A2) = (-l)''-^|J,-M;Ai,A2) . (A3) 

Now rotational invariance of the theory, in the form of 
the Wigner-Eckart theorem, implies that 

{J'M';L'S'\T\J,M;LS) 

= Sj.j,5m,m'{L'S'\T''\LS) 

= (J', -M'; L'S'\T\J, -M; LS) . (A4) 

(J'M';A;,A^|r|J,M;Ai,A2) 

= ^j,j"5m,m'(A'i, AglT'^jAi, A2) 

= (J',-M';A'i,A'2|r|J,-M;Ai,A2) . (A5) 

It follows from Eq. (A3), Eq. (A4) and Eq. (A5) that 
matrix elements between time-reversed states obey 

(ei|r|e/) = mf) , m 

where \i) and |/) are angular momentum eigenstates. It 
then immediately follows from Eq (A2) that the T-matrix 
is symmetric in this basis, 

{f\T\i) = {i\T\f} . (A7) 



[1] E. Komatsu et al. (WMAP), Astrophys. J. Suppl. 192, [2] M. Ackermann et al, JCAP 1005, 025 (2010), 
18 (2011), arXiv: 1001. 4538 [astro-ph.CO]. 



10 



A. Pinzkc, C. Pfrom- 
arXiv: 1105.3240 [astro- 



arXiv: 1002.2239 [astro-ph.CO]; 
mcr, and L. Bcrgstrom, (2011), 
ph.HE]. 

P. Scott et al., JCAP 1001, 031 (2010), arXiv:0909.3300 
[astro-ph.CO]; A. A. Abdo et al., Astrophys. J. 712, 147 
(2010), arXiv:1001.4531 [astro-ph.CO]. 
A. Geringer-Sameth and S. M. Koushiappas, (2011), 
arXiv: 1108.2914 [astro-ph.CO]; M. Ackcrmann et al. 
(Fermi-LAT), (2011), arXiv:1108.3546 [astro-ph.HE]. 
D. Hooper and T. Linden, (2011), arXiv:1110.0006 
[astro-ph.HE]. 

M. R. Buckley and D. Hooper, Phys. Rev. D82, 063501 

(2010), arXiv: 1004. 1644 [hcp-ph]. 

M. CircUi, P. Panci, and P. D. Scrpico, Nucl. Phys. 
B840, 284 (2010), arXiv:0912.0663 [astro-ph.CO]. 
M. Papucci and A. Strumia, JCAP 1003, 014 (2010), 
arXiv:0912.0742 [hep-ph]. 

K. N. Abazajian, P. Agrawal, Z. Chacko, and C. Kilic, 
JCAP 1011, 041 (2010), arXiv:1002.3820 [astro-ph.HE]. 
A. A. Abdo et al. (Fermi-LAT), JCAP 1004, 014 (2010), 
arXiv: 1002.4415 [astro-ph.CO]; G. Hutsi, A. Hektor, 
and M. Raidal, ibid. 1007, 008 (2010), arXiv: 1004.2036 
[astro-ph.HE]; C. Arina and M. H. G. Tytgat, ibid. 1101, 
Oil (2011), arXiv: 1007.2765 [astro-ph.CO]. 
A. A. Abdo et al., Phys. Rev. Lett. 104, 091302 (2010), 
arXiv: 1001. 4836 [astro-ph.HE]. 

L. Bergstrom and P. Ullio, Nucl. Phys. B504, 27 (1997), 

arXiv:hep-ph/9706232. 

Z. Bern, P. Gondolo, and M. Perelstein, Phys. Lett. 

B411, 86 (1997), arXiv:hep-ph/9706538. 

P. Ullio and L. Bergstrom, Phys. Rev. D57, 1962 (1998), 

arXiv:hcp-ph/9707333. 

F. Boudjcma, A. Scmcnov, and D. Tcmcs, Phys. Rev. 
D72, 055024 (2005), arXiv:hep-ph/0507127. 
L. Bergstrom, T. Bringmann, M. Eriksson, and 
M. Gustafsson, JCAP 0504, 004 (2005), arXivihep- 



[17 

[is; 

[19 
[20 

[21 

[22; 

[23; 

[24' 
[25^ 
[26' 
[27 

[28; 
[29; 

[30 
[31 
[32; 



ph/0412001. 

G. Bertonc, C. B. Jackson, G. Shaughnessy, T. M. P. 
Tait, and A. Vallinotto, Phys. Rev. D80, 023512 (2009), 
arXiv:0904.1442 [astro-ph.HE]. 

A. Birkedal, A. Noble, M. Perelstein, and A. Spray, Phys. 

Rev. D74, 035002 (2006), arXiv:hep-ph/0603077. 

M. Perelstein and A. Spray, Phys. Rev. D75, 083519 

(2007), arXiv:hcp-ph/0610357. 

M. Gustafsson, E. Lundstrom, L. Bergstrom, and 
J. Edsjo, Phys. Rev. Lett. 99, 041301 (2007), arXiv:astro- 
ph/0703512. 

C. B. Jackson, G. Servant, G. Shaughnessy, T. M. P. Tait, 
and M. Taoso, JCAP 1004, 004 (2010), arXiv:0912.0004 

[hep-ph]. 

J. Goodman et al, Nucl. Phys. B844, 55 (2011), 
arXiv: 1009.0008 [hep-ph]. 

A. Badin, G. K. Yeghiyan, and A. A. Petrov, (2009), 
arXiv:0909.5219 [hep-ph]. 

M. Jacob and G. C. Wick, Ann. Phys. 7, 404 (1959). 

L. D. Landau, Dokl. Akad. Nawk., USSR 60, 207 (1948). 

C.-N. Yang, Phys. Rev. 77, 242 (1950). 

G. F. Giudicc, R. Rattazzi, and J. D. Wells, Nucl. Phys. 

B544, 3 (1999), arXiv:hep-ph/9811291. 

T. Han, J. D. Lykken, and R.-J. Zhang, Phys.Rev. D59, 

105006 (1999), arXiv:hep-ph/9811350 [hep-ph]. 

R. Catena and P. Ullio, JCAP 1008, 004 (2010), 

arXiv:0907.0018 [astro-ph.CO]. 

A. A. Abdo et al. (The Fermi-LAT), Phys. Rev. Lett. 
104, 101101 (2010), arXiv: 1002.3603 [astro-ph.HE]. 
K. N. Abazajian, S. Blanchet, and J. P. Harding, (2010), 
arXiv:1011.5090 [hep-ph]. 

A. Galper, R. Aptekar, I. Arkhangelskaya, M. Boezio, 
V. Bonvicini, et al, Astrophys. Space Sci. Trans. 7, 75 
(2011). 



11 



